Quantum transport through a deformable molecular transistor 
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The linear transport properties of a model molecular transistor with electron-electron and electron- 
phonon interactions were investigated analytically and numerically. The model takes into account 
phonon modulation of the electronic energy levels and of the tunnelling barrier between the molecule 
and the electrodes. When both effects are present they lead to asymmetries in the dependence of 
the conductance on gate voltage. The Kondo effect is observed in the presence of electron-phonon 
interactions. There are important qualitative differences between the cases of weak and strong 
coupling. In the first case the standard Kondo effect driven by spin fluctuations occurs. In the 
second case, it is driven by charge fluctuations. The Fermi-liquid relation between the spectral 
density of the molecule and its charge is altered by electron-phonon interactions. Remarkably, the 
relation between the zero-temperature conductance and the charge remains unchanged. Therefore, 
there is perfect transmission in all regimes whenever the average number of electrons in the molecule 
is an odd integer. 

PACS numbers: 71.27.+a,75.20.Hr,73.63.-b,85.65.+h 



I. INTRODUCTION 

There has been in recent years a surge of interest in 
the study of transport phenomena in nanoscale systems 
motivated by their potential as electronic devices. These 
systems differ fundamentally from conventional conduc- 
tors in that their electronic spectrum is discrete and their 
charging energies can not be ignored. 

Charging effects on the transport properties have been 
intensively studied. Coulomb blockade has been ob- 
served in the conductance of quantum dotsi and in single- 
molecule devices weakly coupled to the electrodesi^ 
The enhancement of the low-temperature conductance in 
the valley between Coulomb blockade peaks due to the 
Kondo effect has also been observed both in quantum 
dots^ and in single molecules having well defined spin 
and charge statesi 2 i 5 i 6 i 7 

An important feature of molecules is that they gen- 
erally distort upon the addition or the removal of elec- 
trons. While in conventional systems the energies asso- 
ciated with atomic motion are much lower than typical 
electronic energies, this is not necessarily true in molecu- 
lar devices. It has been recently shown that the Coulomb 
charging energies of single molecules can be considerably 
reduced by screening due to the electrodes^ The former 
may drop from a few eV, for an isolated molecule, to a 
few hundreds of meV becoming then of the same order 
of magnitude than vibrational energies. Interesting new 
physics emerges when the two energy scales become com- 
parable as was observed in inelastic electron tunnelling 
spectra of small molecules adsorbed on surfaces^ in Ceo 
molecular-scale transistors^ and in a suspended quantum 
dot cavity^ 

In this paper we study the linear transport properties 
of a model molecular device in which electron-electron 
and electron-phonon interactions are present. Previous 
authors investigated similar models in the regimes of high 
temperature or weak electron-phonon coupling. iLiSii^ 
We are interested here in the low-temperature transport 



properties in regimes for which the charging energy and 
the electron-phonon energies can be of the same order 
and large compared to the broadening of the electronic 
levels. This requires treating both interactions on the 
same footing and non-perturbatively. 

In a vibrating molecule connected to electrodes the po- 
sition of the electronic energy levels with respect to the 
Fermi level and the height of the tunnelling barrier be- 
tween the molecule and the electrodes are phonon modu- 
lated. These effects have important consequences for the 
transport. In a previous publication^ we gave a short 
account of our results for a model in which only the first 
effect was included. Here, we present a detailed analysis 
of the model, including the effects of tunnelling-barrier 
modulation. 

Our results are analytical and numerical. Fermi-liquid 
theory was used to derive some general properties of the 
zero-temperature conductance and detailed information 
on its dependence on the parameters of the model was ob- 
tained in the limits of weak and strong electron-phonon 
coupling. Numerical calculations were also performed in 
all the relevant parameter regimes using the Numerical 
Renormalization Group method at zero and at finite tem- 
perature. 

Our main results are as follows. In all parameter 
regimes there is a peak in the zero-temperature conduc- 
tance as a function of gate voltage G(V g ) whose hight cor- 
responds to perfect transmission through the junction. It 
occurs at the value of V g for which the average number 
of electrons in the device is an odd integer. 

The width and shape of the peak depend on the type 
and strength of the electron-phonon coupling. When ei- 
ther energy level modulation (ELM) or tunnelling barrier 
modulation (TBM) are present (but not both), G(V g ) is 
symmetric around its maximum. However, if the two are 
simultaneously present, the curve is asymmetric. This 
feature is very pronounced in some cases. The two cou- 
plings have opposite effects on the width of G(V g ). ELM 
leads to peak narrowing while TBM has the opposite ef- 



feet. 

There are parameter regimes in which resonant trans- 
mission is due to the Kondo effect. When the ampli- 
tude of ELM ranges from low to moderate, the effect 
is qualitatively similar to that observed in the absence 
of the electron-phonon coupling. When ELM is strong, 
however, the nature of the ground state of the molecule 
changes and the Kondo effect does not result from spin 
fluctuations (as is the case in quantum dots) but from 
charge fluctuations. There are important qualitative dif- 
ferences between the two cases. 

At finite temperature, in the Kondo regime, the height 
of the conductance peak decreases and is strongly sup- 
pressed beyond a characteristic temperature, the Kondo 
temperature. An increase in TBM by itself leads to an 
increase of the Kondo temperature. The effect of in- 
creasing ELM is non monotonic: the Kondo temperature 
increases for small coupling but it decreases for strong 
coupling. 

Above the Kondo temperature well defined asymmetric 
Coulomb blockade peaks are observed if ELM is not too 
strong. When its strength increases, the Coulomb peaks 
become closer to each other and disappear in the strong 
coupling regime. The effect of TBM is to smear these 
peaks when they exist. 

The rest of this paper is organized as follows. In Sec- 
tion [H] we describe the model that we use in our calcula- 
tions. Section IlIII contains our analytical results, includ- 
ing the derivation of the Fermi-liquid relations for our 
model and the analysis of the limiting cases of weak and 
strong electron-phonon coupling. Our numerical results 
are presented and discussed in Section IIVI Finally, we 
state the conclusions of our study in Section Ivl 



II. THE MODEL 

We model a molecule connected to metallic electrodes 
in a range of gate voltages such that only one molecu- 
lar orbital effectively participates in the transport. The 
molecule has a symmetric vibrational mode of frequency 
ujq coupled to the electronic coordinates. The energy 
e d {x) of the molecular level and the tunnelling matrix el- 
ements between the molecule and left (L) and right (R) 
electrodes, Vg(x) = V(x) {I — L,R), depend on a dimen- 
sionless vibrational coordinate x. For small distortions 
these quantities may be expanded as e d (x) « e d — A x 
and V(x) w V [1 + g x] where A and g are two coupling 
constants^ The first one is an energy scale, the second 
one is dimensionless. A scheme of the device is shown in 
Fig. © 

The Hamiltonian of the combined system consists of 
three terms that describe the isolated molecule, the iso- 
lated electrodes, and the coupling between them, respec- 
tively, 

H = Hm + He + Hm-e ■ (1) 
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FIG. 1: A scheme of the model device studied in this paper. A 
molecule of average length Lo is connected to source and drain 
electrodes. As it vibrates symmetrically, the distance between 
the end groups and the electrodes fluctuates thus modulating 
the tunnelling barrier between the elements. There is also a 
modulation of the position of the molecular energy levels with 
respect to the Fermi level not represented in the scheme. 

The Hamiltonian of the molecule is 

H M = e d h d + U n d pi di - A (n d - I) (a + a 1 ) + oj a^a ,(2) 

where h d = d\d a is the charge operator and U is the 
Coulomb repulsion between two electrons on the same 
molecular orbital. In appropriate units we write the elon- 
gation of the molecule as x = a + , where a and are 
the phonon operators. 

The Hamiltonian of the leads is 

He = J2 £ ^ c Li c ^ = L, R) , (3) 

where we consider for simplicity the case of identical elec- 
trodes with dispersion £R(k) = ei/k) = e(k). We assume 
for convenience that the conduction band is symmetric 
with respect to the chemical potential and we take the 
latter as the origin of energies. 

The third term in Eq. is the tunnelling Hamiltonian 
between the molecule and the electrodes, 

H M -e = [l+g{a + <J)]Vj2 (4<W + h.c) (4) 

= [l + gia + a^Vv , 

where we assumed that the tunnelling matrix element 
V is k-independent and we defined the operator v for 
notational convenience. 

Note that if either g — or A = and e d = —U /2, the 
Hamiltonian of the model is symmetric under an electron- 
hole transformation (plus the inversion x — ► —x if A 0). 
At the symmetric point, the average number of electrons 
of the molecule is n d = (n d ) — I. If both g and A are 
non vanishing, this symmetry is lost. We shall see below 
that this has important consequences. 

An expression for the conductance G of the molecu- 
lar junction at zero bias can be derived using standard 
methodsi 16 ' 17 We find 

(5) 
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where f(u>) is the Fermi function, A = 2ttpoV 2 and po 
is the electronic density of states of the electrodes eval- 
uated at the Fermi level. The spectral density f>dd{oj) = 
— n ImGdd(^), where we defined a modified d-electron 
Green function Gdd{^>) as 



G dd (oj) = -i dte l ^([{l + gx{t))dv{t), 



(1 + gx) 4] + ) 



(6) 



III. ANALYTICAL RESULTS 
A. Exact results at zero temperature 

The ground state of Hamiltonian JQ) is expected to be a 
Fermi liquid. Exact relationships can then be established 
between the zero-temperature conductance through the 
junction, the spectral density of the molecule at the Fermi 
level and its charge. To derive these relations the struc- 
ture of the expansion of the electronic Green functions 
in powers of the electron-electron and electron-phonon 
interactions coupling constants must be examined. 

We start by noting that when the molecule is coupled 
to the contacts it will in general deform. It is convenient 
to carry out the perturbation expansion around the de- 
formed ground state. We thus apply a canonical trans- 
formation that shifts the phonon operators, a — > b + a, 
and write 



#(a,a + ) = H{b + a,tf +a) = H' a {b,tf) 



which defines the transformed Hamiltonian H' a (b, b'). 

Since the ground state energy is invariant under this 
transformation, 



E ee (H) = (H' a ) , 



(8) 



we may take the derivative with respect to a of this equa- 
tion to obtain, 

lo (b + & f ) Q + 2aw - 2A(n d - l) a + 2gV(v) a = , (9) 

where (• • -) Q denotes an expectation value with respect to 
the exact ground state of the transformed Hamiltonian. 
The deformation (x) of the molecule is found by imposing 
that (b + b^a — on the transformed ground state: 

(x) A , , q , , . . 

a= V = — ("d-l)a- — V(v) a . (10) 

The transformed Hamiltonian can be expressed as a sum 
of one-body and many-body terms as 



H — Hq + H [ , 



H' a = H E +u tfb + e d h d + V v 



(11) 



(12) 



H[ = (b + tf)[gV(v 
+ Un d ]Udi , 



(v) a ) - A (fl d - (hd}a)} (13) 



where we dropped a constant energy shift and we defined 
the renormalized parameters 



£d = £d — 2Aa , and V = V (1 + 2ga) 



(14) 



The molecule is only coupled to a local symmetric lin- 
ear combination of the states of the leads, 



|tf c ) =N- l/2 ' 



I* 



k.R/ 



k.L/ 



k 



V2 



(15) 



where TV is a normalization factor. The remaining 
conduction-electron states can be integrated out reduc- 
ing the problem to an effective interacting two-site model 
with a damping term that reflects the coupling of the 
subsystem to the electrodes. Physical quantities may be 
obtained from the 2x2 Green-function matrix 



G cc (lo) Gcd(u) 
G dc {u) G d d{u) 



(16) 



The Green function associated to the one-body Hamilto- 



nian H' a is 



GoM = 



1/G° C H -V2V 



-V2V 



(17) 



w - (d 



where, in the wideband limit, G° c (w) = J2u l/( tJ ~ £ k) = 
—iirpQ and po is the bulk density of states of the elec- 
trodes. Note that Gq(uj) is not the non-interacting Green 
function and it includes effects of the interactions to all 
orders through the renormalized parameters (|14|) . 

The remaining effects of the interaction H[ are em- 
bodied in a self-energy matrix defined through the Dyson 
equation 



G" 1 ^) = GoV) - S(w) 



(18) 



In a Fermi-liquid ground state all the elements of the 
self-energy matrix are purely real at the Fermi level. We 
checked that this is the case for our model by computing 
a few low-order diagrams in the expansion of S in powers 
of H[. Luttinger's theoremiSiiScan then be generalized 
to the present situation in order to establish some useful 
Fermi-liquid relationships. 

We found a simple relationship between p cc (0), the c- 
component of the spectral density at the Fermi level, and 
the electronic population of the molecule in the physically 
relevant wideband limit. It reads, 



p cc (0) = pocos 2 {^ndj 



(19) 



Equation l|19l) is an exact result valid for all values of 
the parameters of the Hamiltonian (JIJ. For 5 = but 
keeping A arbitrary we found an equivalent expression 
for p d d{0), the projection of the spectral density on the 
molecular orbital. In the wideband limit this is 



Pdd(0)\ c 



1 

7rA 



»d 



(20) 
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The content of Eq. I|20[l is that, for g = 0, the d-spectral 
density of the interacting system at the Fermi level is 
pinned at the value it takes for a non-interacting system 
with the same electron occupancy. This is a well known 
result for model (JTJ in the absence of electron-phonon 
coupling^ We see that it remains valid for all values of 
A and g — 0. 

However, Luttinger's theorem does not lead to a sim- 
ilar result for pdd hi the general g ^ case and the ele- 
ments of the self-energy matrix appear explicitly in the 
corresponding expression. 

There is some simplification in the particular case A = 
0, ed — —U/2. Then, the Hamiltonian (JTJ is electron 
hole symmetric [c.f. Section [H], na = 1, and T, cc (0) and 
Sd<j(0) vanish. In this case we find 



B. Analysis of limiting cases 

Many features of the solution of our problem can be 
found from an analysis of the limiting cases of small and 
large electron-phonon coupling. 

We start from the Hamiltonian of the isolated molecule 
given in Eq. (J2J). Its eigenvalues and eigenfunctions can 
be written down explicitly^ 

\0,m) = U+\0) e \m), £ ,m = ~ + mcj , 

\<r,m) = \a) e \m) , E a . m = e d + muj , (27) 

|2, m) = U~ |U) e \m) , E^ rn = + 2e d + U + ™„, 



Pdd(Q)\ X =o,n d =l 



7rA 



(21) 



where the subscript e denotes electronic states, |m) is 
TO-th excited state of the harmonic oscillator and 



where 



(22) 



It can be easily shown that S c y is of order g 2 and higher. 
It then follows from Eqs. l(Ttj)). |Q3J) and that A = 
A (l + 0(.g 2 ) + ■•■). 

Interestingly, the non-universality of pdd(0) (of which 
Eq. 12111 is only an example) does not extend to the T = 
conductance. The latter is obtained setting T — in 
Eq. ©. 



= ttA p dd (0) 

LrQ 



(23) 



where Go = 2e 2 /h is the quantum of conductance. 

An expression for the modified Green function Gdd{u) 
can be obtained by writing down the equation of motion 
for the conduction-electron Green function G cc (lo) using 
the definition © and the original Hamiltonian QJ. We 
find 

G CC ( W ) = G2» + 2V 2 G° C HG M MG° C M . (24) 



U ± = exp 



±— (a - at) 

LUq 



(28) 



At Ed = —U/2 the states with zero and two electrons 
are degenerate. We take this as a reference point and 
write 



U 

^ + V d 



(29) 



The number of electrons in the ground state of the 
molecule is either odd or even depending on whether 
2A 2 /wo is smaller or larger than U. It is convenient to 
define an effective interaction parameter, 



2A 2 
w 



(30) 



It will be seen below that the physics in the cases U e g > 
and f/gg < is quite different. 



1. 



U eff >0 



Setting lu — and taking the imaginary part of the above 
equation this becomes 



Pcc(0) = pa [1 - 7rAp d(i (0)] . 



(25) 



Comparing Eq. (|25() with the Fermi-liquid relationship 
(fT§|) we find 



TrApddiO) = q^ = Slir {2 



(26) 



The zero-temperature conductance thus depends only on 
the occupation of the molecular orbital just as in the 
absence of tunnelling barrier modulation and it reaches 
its maximum value Go when rid — I. A related result 
was obtained in Ref. [ll] in perturbation theory in A for 
U = and g = 0. 



In this case the ground-state of the isolated molecule 
for I Vd\ < U/2 is the spin-doublet \a, 0) . There is a charge 
excitation gap given by ~ U c s = U — 2X 2 /u>o for Vd = 
and a phonon gap ojq. The low-energy excitations are 
spin fluctuations. 

An effective Hamiltonian for excitations involving the 
spin doublet can be derived by perturbatively project- 
ing out the empty and doubly occupied states from the 
space of available states by means of a Schrieffer- Wolff 
transformation. To lowest order in V/U, the low-energy 
effective Hamiltonian for spin fluctuations H e ff is 



-ffeff = ujQa^a + He + JS • s c 



(31) 



where S is the spin operator for the doublet and s c = 
1/2 , c| <t s ,s' <V is the spin operator of an electron 
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on the state (|15fl coupled to the molecule. The coupling 
constant J is given by 



2. 



J=^f] V K0|[l + .g(« + at)]^|m)| 2 

^ m— }J-= ± 



2A 2 /w ^ + 2mw /£7 



(32) 



where Jo = 8A/(7r£/po) and we have taken Vd — for 
simplicity. 

A simple analytical expression can be derived in the 
limit of small electron-phonon coupling, g, X/U <C 1. We 
find2a 



J = Ja 



1 



(2X/U) 2 +g 2 
1 + 2lu /U 



(33) 



Therefore, J increases with the electron-phonon coupling 
in the weak coupling limit. 

Equation (|31() is the well known Kondo model Hamil- 
tonian , 19 The dependence of the conductance on temper- 
ature, gate- voltage and magnetic field is well understood 
in the absence of electron phonon coupling. We summa- 
rize below the main features for g = A = 0. 

At T = 0, there is a narrow resonance of width Tk ~ 
Dexp (— l/(Jp )) in the d-electron spectral density at 
the Fermi level. This resonance provides a channel for 
conduction and, at the symmetric point, rid — 1, G = Gq 
in agreement with Luttinger's theorem. 

For Vd 7^ 0, the gap between the ground state doublet 
and the empty state (Vd > 0), or the doubly occupied 
state (Vd < 0) is reduced. The charge then deviates from 
rid = 1< The Kondo resonance shifts with respect to the 
Fermi level and the pdd(0) decreases. For \Vd\ > U/2 
the Kondo resonance disappears. There is thus a peak in 
G(V g ) whose width ~ U is large. An applied magnetic 
field B also destroys the Kondo resonance by breaking 
the symmetry between the spin up and spin down ground 
states. This happens for B ~ Tk which is very small in 
the Kondo limit. Therefore, the peak in G(B) is narrow. 
At finite temperature the conductance at the symmetric 
point decreases becoming very small for T > Tr. In this 
temperature range G exhibits Coulomb blockade peaks 
at gate voltages Vd ~ ±U/2. 

These features are present throughout the weak cou- 
pling regime with some minor modifications. It follows 
from equation (|33|l that Tk increases with increasing 
A and g. Therefore, the peak in G(B) becomes wider 
and the temperature variation of the conductance at 
resonance becomes slower. If both g and A ^ the 
zero-temperature peak in G(Vd) shifts from Vd = to 
Vd = V^ such that nd = 1 and acquires an asymmetric 
shape. These features are related to the loss of particle- 
hole symmetry referred to above. The conductance is 
higher to the left of the maximum than to its right be- 
cause for negative values of Vd — VJ the doubly occupied 
state that has a stronger coupling to the electrodes is fa- 
vored. The width of the peak and the separation between 
the finite-temperature Coulomb blockade peaks, that are 
also asymmetric, decrease with increasing A mirroring the 
evolution of U„g. 



U eff <0 



In this regime the ground-state for Vd = is a charge 
doublet formed by the states |2, 0) and |0, 0). The low- 
energy excitations of the system are thus charge fluctu- 
ations and there is a large gap ~ 2A 2 /wo — U for spin 
excitations. 

We discuss the case g = first. A low-energy ef- 
fective Hamiltonian can be found as before applying a 
Schrieffer- Wolff transformation, but now we project out 
the singly occupied states. We introduce pseudo-spin op- 
erators and r* with eigenvalues ±1/2, corresponding 
to the eigenvalues Q = 2 and Q = of the charge of the 
molecule and of the orbital c, respectively. We also define 

raising and lowering operators , = (tJ~) = cijcij and 
similar ones for the leads. In terms of these, the effective 
Hamiltonian for Vd = g = is: 



H = H E + J,, rj tI 
with the couplings 



J± 



T rf +T c + T d T c + ) . 



(34) 



J\\ = Jo E 



"p\U+\m)(m\U-\0) 
' Q 2A 2 /w t/~ 1 + 2muJo/U 



, (35) 



^ (0\U+\m)(m\U+\0) 
- J »L 2A 2 M,t/ - 1 + 2mu; /U " W 

Equations H34fl - H36fl define the Hamiltonian of an 
anisotropic Kondo model (AKM). 20,2 ii The couplings l|35|) 
and JSSJl can be easily estimated in the limit X/loq 3> 1 
by noting that 

(0\U+\m)(m\U-\0) = e -(V-o) 2 (*M)) 2m (3?) 

m! 

= (-l) m (0|[/ + |TO)(m|C/ + |0) . 

Considered as a function of m, the last expression on the 
right hand side of the upper line of Eq. 137|l is strongly 
peaked at m* — (X/uoq) 2 ^> 1, while the denominators in 
Eqs l|35(l and l|36|l are slowly varying functions of m. In 
the strong coupling limit these can be approximated by 
their value at m* and taken out of the sum. This results 



Jii 



A 2 



Ji 



J\\e 



-2(A/u;o) 2 



(38) 



This strong anisotropy originates in the fact that the 
phonon ground states corresponding to the two charge 
states Qd — 2 and Qd — are very different. In the 
strong-coupling limit, the anisotropy ratio J±/J\\ is pre- 
cisely the overlap (2, OjO, 0) between them [c.f. Eq. J2ZJ]. 
The Kondo temperature of the AKM ia 20 i 21 



Takm 



J± 
Jii 



J||P0 



D exp 



TTU>0 ( X 



, (39) 
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where the last expression is its asymptotic form for large 

In contrast with the weak-coupling situation, Takm 
decreases sharply when A increases. Application of a 
small potential \Vd\ > ?akm splits the charge doublet 
and destroys the Kondo resonance. Hence, the peak 
in G(Vd) is very narrow in this limit. Conversely, the 
doublet is insensitive to a magnetic field and the width 
~ 2A 2 /cjo — Uot the peak in G(B) is now large. Finally, 
no Coulomb blockade peaks are to be seen upon applica- 
tion of a gate voltage because it leads to an increase of 
the energy difference between a ground state with zero 
or two electrons and an excited state with one electron 
irrespective of its sign. These features are opposite of 
those that characterize the properties of G in the weak 
coupling regime. 

We discuss now the case g ^ 0. The modulation of the 
tunnelling amplitude has two effects. The first one is to 
renormalize the coupling constants, 



J,, -> Jjj = J,, [1 + (gX/tuo) 2 ] , 



(40) 



J± = J±[i- {gX/u f 



The second effect is the appearance of a term in the 
Hamiltonian that breaks the symmetry between the two 
charge states, favoring the doubly occupied one: 



^AKM — > ^AKM — ~ J\\ T - 
UJ 



d i 



(41) 



where -ff^KM denotes the AKM Hamiltonian with the 
new couplings. These effects can be understood quite 
easily. In the state with charge Qd — 2 the molecule elon- 
gates. This increases its coupling to the electrodes and 
its energy is lowered with respect to that of the empty 
state. The tunnelling barrier between them also increases 
which translates into a reduction of Jj_ . 

We can use the modified couplings (j40|) to compute the 
(/-dependent Kondo temperature. The result is, 



rpl 

AKM 

?akm 



exp 




(42) 



In the strong coupling regime the Kondo temperature in- 
creases or decreases with g, depending on the ratio X/luq. 

The consequences of a finite g on the conductance are 
twofold. First, the peak in G(Vd) no longer occurs at 
Vd = but it shifts to ~ g Jn X/wq, the value for 
which the asymmetry generated by the second term in 
Eq. (|4*T1) is compensated. Secondly, as in the previous 
case, the singly occupied states are more strongly coupled 
to the doubly occupied state leading to asymmetry in the 
shape of G(Vd). This time, the asymmetry is much more 
pronounced than for weak coupling because the relative 
weights of the components of with Qd = and Qd = 2 
in the wavefunction are more sensitive to Vd- 



IV. NUMERICAL RESULTS 

In this Section we present the numerical solution of 
model Q using the numerical renormalization group 
method (NRG)—£&2^2^ incorporating a modification de- 
signed to improve the accuracy of the computation of the 
spectral density^ 

Results for the case g = were discussed in a previous 
paper^* Here, we focus on the effects of tunnelling barrier 
modulation on the conductance. 

In our numerical calculations we used half the band- 
width of the conduction electron band D as the unit of en- 
ergy. The parameters U — 0.1, A = 0.01, and luq = 0.05 
were kept fixed and we studied the conductance as a func- 
tion of g, X and Vd- 



A. The electron hole symmetric case at A = 

We start by analyzing the electron-hole symmetric 
case, A = 0, V d = [c.f. Eq. Pty]. Fig. El shows the 
spectral density pdd{u) at zero temperature for several 
values of g between and 0.5. 

The top curve corresponds to g = 0. It exhibits a 
narrow Kondo resonance centered at the Fermi level and 
two Coulomb peaks at ui ~ ±U/2. The height of the 
Kondo resonance is 1/irA ~ 32 in this case. 

With increasing g the hight of the central peak de- 
creases and its width increases. This evolution is shown 
in more detail in the left inset to the figure. 

That the height of the central peak depends upon g 
shows that pdd(0) does not obey Luttinger's theorem in 
the simple form of Eq. H20|) . An effective hybridization 
parameter A may be obtained from the numerical values 
of Pdd{0) using Eq. J2J. Its dependence on g is repre- 
sented in the right inset to the figure. It is seen that A 
increases quadratically with g in the electron-hole sym- 
metric case as was anticipated in Section llll Al 

The figure shows that the width T of the Kondo res- 
onance increases monotonically with g. The dependence 
of r on g is represented in the inset to Fig. [5] that shows 
that logT varies quadratically with g for small g. This 
behavior can be understood from Eq. that predicts 
a quadratic increase of the Kondo coupling J for small 
g and A = 0, recalling that T ~ Tk and that the latter 
varies exponentially with J . 

Finally, we observe that the Coulomb peaks broaden 
with increasing g until they eventually disappear when 
the system enters the mixed- valence regime for large g. 

The temperature dependence of the conductance at 
Vd = is shown in Fig. [31 for the same values of g as in 
Fig. [21 G(g, T) reaches its maximum value Gq a,t T = 
for all g in agreement with Eq. H2(j|) ■ This confirms the 
result of Fermi-liquid theory that tunnelling barrier mod- 
ulation does not lead to a renormalization of the value of 
the conductance at resonance despite that the spectral 
density at the Fermi level is modified. 
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FIG. 2: Spectral density p dd (ui) at zero temperature in the 
electron-hole symmetric case. The electron-phonon coupling 
g varies between and 0.5 in steps of 0.05 from top to bottom 
in the region of the central peak. Left inset: a zoom of the 
central part of the main figure. Right inset: g-dependence of 
the effective hybridization parameter A. 



FIG. 4: Spectral density p dd (uj) at T = for A = -0.06 and 
g = 0.2. Thick line: V d /U = V d * /U = 0.159, n d = 1. Dashed 
line: V d /U = 0.163, n d < 1. Thin line: V d /U = 0.155, n d > 1. 
Inset: a zoom of the central peak showing the evolution of the 
Kondo peak with varying V d . 




T/D 

FIG. 3: Temperature dependence of the conductance in the 
electron-hole symmetric case. The electron phonon couplings 
are as in Fig. [2] The curves correspond to increasing values 
of g from left to right. Inset: g-dependence of the width of 
the central peak of the curves shown in the previous figure 

The conductance at resonance decreases with tempera- 
ture. The temperature scale for this process is the Kondo 
temperature. The observed shift of the curves to higher 
temperatures with increasing g reflects the evolution of 
the latter. 



B. The general case, A and g ^ 

We pointed out in Section [H] that when both A and g 
are non zero the system looses electron-hole symmetry. 
This is at the origin of several observable features in the 
frequency dependence of the molecule's spectral density 



and in the dependence of the conductance on gate volt- 
age. All the results described below were obtained with 
.9 = 0.2. 

Figure 0] shows the (i-electron spectral density at zero 
temperature in the strong coupling regime, U q q < 0. We 
display results for three values of V d : V£, chosen so that 
n d = 1, and V d = ± SV d . 

The thick line represents pdd{u) for V d — V d - We 
observe a very narrow peak at the Fermi level due to 
the anisotropic Kondo effect described in Section lill B 21 
We also observe sideband peaks that are associated to 
excitation of the higher levels of the isolated molecule 

at uj ~ i^eff [ c '^' ^ ^ s seen these are 

not positioned symmetrically with respect to the central 
peak and that their widths and heights are different. As 
discussed above, the origin of the asymmetry is that the 
states with Q d = and Q d = 2 are coupled with different 
strength to the states with Q d = 1. The matrix element 
for the transition |2, 0) — > |1<7, 0) is larger than that for 
the transition |0, 0) — * |1<7, 0). 

The spectral densities for V d — V d ± SV d are repre- 
sented in the figure by dashed and solid thin lines, re- 
spectively. For our choice of parameters, \6V d \ < ?akm- 
Therefore, the Kondo resonance continues to exist but it 
shifts and its amplitude decreases as shown in the inset 
to Fig. 0| which is a detailed view of the central part of 
the main plot. 

When V d > V rf * the empty orbital has a larger weight in 
the ground state wavefunction than the doubly occupied 
orbital. Then, the sideband at u> > is enhanced and 
that at u> > is suppressed. The opposite occurs for 

Vd < vi 

For \5V d \ > Takm ( no * snown )j the molecule becomes 
fully charge-polarized^ and the Kondo peak is com- 
pletely suppressed. 

We have also computed the spectral density in the 
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FIG. 5: Conductance as a function of Vd and temperature 
for weak coupling, J7 e gp > 0. Panel (a): A = 0.02. Panel 
(b): A = 0.035. Temperatures are: T = 0,7.3 x 10~ 4 , 1.5 x 
10" 3 ,3 x 10" 3 ,6 x 10 -3 , 0.023 and 0.094 from top to bottom. 
In all cases g = 0.2. 



weak coupling regime, U e Q > (not shown) . For moder- 
ate values of g, Vj ~ and the effects of the asymmetry 
are much less pronounced than for the strong coupling 
case just discussed. 

We now turn to the analysis of the results for the con- 
ductance. Fig. [S] shows G(Vd) for several temperatures 
and two values of A, both in the regime U e Q > 0. We ob- 
serve a broad conductance peak of width ~ f7 e £p at zero 
temperature. There is perfect conductance at V? ^ 
and the peak is asymmetric with respect to Vi. 

At finite temperature the Kondo effect is first sup- 
pressed at the center of the peak where the Kondo tem- 
perature is the lowest, leading to a two-peak structure. 
There are two Coulomb blockade peaks whose widths and 
hights are different for the same reasons that make the 
spectral density asymmetric. 

A comparison between Figs. EJa) and [5fb) shows 
that an increase of A leads to an enhancement of the 
anisotropy. This is easily understood as we saw that 




-1 -0.5 0.5 1 



V./U 

d 

FIG. 6: Panel (a): Conductance as a function of Vd and 
temperature for {7 e fj = (A = 0.05). Temperatures are 
T = 0, 3 x 10~ 3 , 6 x 10" 3 , 1.2 x 10" 2 , 2.4 x 10 -2 and 4.7 x 10" 2 . 
Inset: Check of the validity of Luttinger's theorem, Eq. 1261 . 
Panel (b) Same as (a) for U e g- < (A = 0.06). Temperatures 
are T=0 and 0.012. Inset. Zoom of the conductance peak for 
temperatures T = 1. 1 x 10~ 5 , 2.3 x 10~ 5 , 4.6 x 10~ 5 , 9.2 x 10~ 5 
and 1.83 x 10~ 4 from top to bottom. For all the curves shown 
in this figure g = 0.2. 

the latter arises from the simultaneous presence of ELM 
and TBM. We also observe a narrowing of the Coulomb 
blockade valley that originates from the decrease of U e Q. 

A further increase of A takes the molecule away from 
the standard Coulomb blockade regime. Fig. shows 
results for U e Q = 0. Here, the differences between the 
energies of the four charge states of the isolate molecule 
are comparable with their width and no Kondo effect is 
expected to occur. No Coulomb blockade peaks appear 
upon raising the temperature as we argued in Section 
IITTTJ^I 

The conductance in the strong-coupling case U b q < 
is shown in Fig. EJb). As in the previous regimes the 
zero-temperature conductance is enhanced and reaches 
the quantum of conductance at some value of Vd but, 
now, the peak is exceedingly narrow. This is the regime 
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where the anisotropic Kondo effect occurs. The width of 
the peak in G(Vd) is now Takm not |£/ e ff|- A very small 
shift \SVd\ > ?akm is sufficient to destroy the Kondo 
resonance. The asymmetry generated by the tunnelling 
barrier is in this case extremely pronounced because the 
charge on the molecule switches from rid = and rid = 2 
within a very small interval of values of Vd- Similarly, 
the conductance at resonance decreases very rapidly with 
increasing temperature as shown in the inset to Fig-EJb) . 

We close this Section with a comparison between the 
conductance of the device calculated from Eq. (JSJ and 
that computed from Eq. 12tj|) . using rid(Vd) obtained from 
the NRG calculation. The data are shown in the inset to 
Fig-EI a )- It is seen that the Fermi- liquid relation l(2*Hjl is 
satisfied within the numerical error. 



V. CONCLUSIONS 

We studied the linear transport properties of a model 
molecular transistor in the presence of electron-electron 
and electron-phonon interactions analytically and by nu- 
merically solving the model using the Numerical Renor- 
malization Group method. 

Electron-phonon interactions lead to modulation of the 
positions of the energy levels of the molecule with respect 
to the Fermi level and of the height of the tunnelling 
barrier between the molecule and the electrodes. These 
effects give rise to new features in the spectral and trans- 



port properties of the system. 

We found that when there is tunnelling barrier modu- 
lation, the strength of the coupling between the molecule 
and the leads depends on the charge of the former. In 
general this effect breaks electron-hole symmetry leading 
to an asymmetric curve of conductance as a function of 
gate voltage. 

The Kondo effect occurs in the presence of electron- 
phonon interactions but its nature is quite different for 
weak and strong coupling. In the first regime the ground 
state of the isolated molecule is a spin doublet and this 
leads to a standard Kondo effect with renormalized pa- 
rameters. In the second regime the lowest lying states of 
the isolated molecule form a charge doublet and the asso- 
ciated charge fluctuations are described by an anisotropic 
Kondo model in a narrow interval of gate voltages. The 
two charge states of the doublet are coupled differently 
to excited states which gives rise to pronounced asym- 
metries in the conductance and in the spectral density. 

We established Fermi-liquid relationships for the inter- 
acting model. Tunnelling barrier modulation leads to a 
non-universal relation between the spectral density of the 
molecule and the electron occupancy. Quite remarkably, 
the relation between the zero-temperature conductance 
and the charge remains unchanged. An important conse- 
quence is that there is perfect transmission in all regimes 
when the number of electrons in the molecule is an odd 
integer. 
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